First-principles studies of the effects of impurities on the ionic and electronic 

conduction in LiFeP0 4 
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Olivine-type LiFePCk is widely considered as a candidate for Li-ion battery electrodes, yet its 
applicability in the pristine state is limited due to poor ionic and electronic conduction. Doping can 
be employed to enhance the material's electrical conductivity. However, this should be understood as 
incorporating electrically active impurities to manipulate the concentration of native point defects 
such as lithium vacancies and small hole polarons which are responsible for ionic and electronic 
conduction, respectively, and not as generating band-like carriers. Possible effects of monovalent (Na, 
K, Cu, and Ag), divalent (Mg and Zn), trivalent (Al), tetravalent (Zr, C, and Si), and pentavalent 
(V and Nb) impurities on the ionic and electronic conductivities of LiFePC>4 are analyzed based on 
results from first-principles density-functional theory calculations. We identify impurities that are 
effective (or ineffective) at enhancing the concentration of lithium vacancies or small hole polarons. 
Based on our studies, we discuss specific strategies for enhancing the electrical conductivity in 
LiFePCU and provide suggestions for further experimental studies. 
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I. INTRODUCTION 

Olivine-type LiFePC>4 has been proposed as a candi- 
date for rechargeable Li-ion battery electrodes because 
of its structural and chemical stabilities, high intercala- 
tion voltage, high theoretical discharge capacity, environ- 
mental friendliness, and potentially low costs^— How- 
ever, poor ionic and electronic conduction are major chal- 
lenges. In addition to nanostructuring and carbon coat- 
ing, doping has been considered as an important path to- 
ward enhancing the electrical conductivity. Since it was 
reported by Chung et al^ that doping LiFePC^ with im- 
purities such as Mg, Ti, Zr, and Nb significantly enhances 
the conductivity, there have been numerous experimental 
works carried out along this path. Nonetheless the role 
of dopants in the conductivity enhancement is still under 
debated— Thus it is important to understand the mecha- 
nisms for ionic and electronic conduction and the effects 
of impurities on the ionic and electronic conductivities. 
Such an understanding is essential in formulating strate- 
gies for improving the electrical (i.e., ionic and/or elec- 
tronic) conductivity and hence the electrochemical per- 
formance of the material. 

Recently, we have carried out comprehensive first- 
principles density-functional theory (DFT) studies of na- 
tive point defects and defect complexes in LiFePC^. 7 
Based on a detailed analysis of the structure, energet- 
ics, and migration of the native defects, we arrived at 
the following main conclusions: (i) Native point defects 
such as small hole polarons {p + ), negatively charged 
lithium vacancies (V^J), negatively charged lithium an- 
tisites (Lip c ), and positively charged iron antisites (Fe^) 
have low formation energies and hence are expected to be 
present in the material. The mobility of Fe^; is low com- 
pared to , supporting earlier suggestions that Fe^; im- 
pedes Li diffusion and hence reduces the electrochemical 



activity!&2 The relative concentrations of these defects 
are, however, sensitive to the experimental conditions 
during synthesis. This suggests that one can suppress 
or enhance certain native defects via tuning the synthe- 
sis conditions, (ii) Native defects in LiFePC^ cannot act 
as sources of band-like electrons and holes, and the ma- 
terial cannot be doped n- or p-type. Any attempt to de- 
liberately shift the Fermi level to the valence-band maxi- 
mum (VBM) or conduction-band minimum (CBM), e.g., 
via doping with acceptors or donors, will result in spon- 
taneous formation of compensating native defects that 
counteract the effects of doping, (iii) The ionic conduc- 
tion occurs via diffusion of whereas the electronic 
conduction proceeds via hopping of p + , which confirms 
earlier suggestions of a polaronic mechanism in the elec- 
tronic conduct ion The ionic conduction is effectively 
one-dimensional along the 6-axis, assuming there are no 
other native defects or extrinsic impurities with low mo- 
bility that block the Li channels; whereas the electronic 
conduction is effectively two-dimensional in the b-c plane. 

It is evident from the conduction mechanisms of 
LiFePC>4 that the ionic and electronic conductivities are 
dependent on mobility and concentration, and hence on 
the migration barriers and formation energies, of and 
p + . This opens the door to manipulating the ionic and 
electronic conductivities of the material via manipulating 
defect concentrations. In fact, LiFePC^ samples with 
high concentrations of and/or p + (and a low con- 
centration of Fe^;) can be obtained in experiment under 
appropriate synthesis conditions. 7 Also, despite the fact 
that LiFeP04 cannot be doped n- or p-type, one can still 
incorporate into the material electrically active impuri- 
ties that are effective in shifting the position of the Fermi 
level away from that determined by the intrinsic native 
defects, slightly toward the CBM (for impurities that ex- 
hibit a donor-like effect) or VBM (for impurities that ex- 
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FIG. 1: Formation energies of intrinsic/native defects V£[ and 
p + (solid blue lines) and of a donor-like impurity (solid red 
line), plotted as a function of Fermi level with respect to the 
VBM. In the absence of extrinsic electrically active impurities, 
the Fermi level of LiFePCU is at ^t" rt where the formation en- 
ergies and thus concentrations of V£[ and p + are equal. When 
the donor-like impurity is incorporated into the material with 
a concentration higher than that of the positively charged 
native defects, the Fermi level is shifted toward the CBM, 
lowering (increasing) the formation energy of (p + ). As 
a result, the ionic (electronic) conductivity will be enhanced 
(reduced). 



hibit an acceptor-like effect). This results in lowering the 
formation energy and hence enhancing the concentration 
of V-T, or p + . Thus "doping" in LiFeP04 should be un- 
derstood as a way to manipulate the concentration of the 
native defects that are responsible for ionic and electronic 
conduction, rather than as generating band-like carriers. 

The mechanism for manipulating the concentrations 
of p + and V{t by incorporating a donor-like impurity is 
illustrated in Fig. [TJ Note that, in the absence of ex- 
trinsic electrically active impurities, the Fermi level of 
insulating materials such as L1FCPO4 is at /x e =A tl e nt , de- 
termined by a charge neutrality condition that involves 
all possible intrinsic/native defects. This condition re- 
quires that defects with different charge states coexist in 
a proportion that maintains overall charge neutrality in 
the material.— ~— Because of the exponential dependence 
of the concentration on defect formation energy (see the 
next section), /x™* is predominantly determined by the 
positively and negatively charged defects with the low- 
est formation energies (In the illustrated example, these 
defects are and p + ). If electrically active impuri- 
ties are incorporated into LiFeP0 4 with concentrations 
higher than that of the native defects, the Fermi level 
may be shifted away from /x™' as the charge neutrality 
condition is re-established) 14 ' 15 as illustrated in Fig. [TJ 

Clearly, it is important to know which impurities are 
effective at enhancing which native defect (s). In this 



Article, we report our first-principles studies of various 
monovalent (Na, K, Cu, and Ag), divalent (Mg and Zn), 
trivalent (Al), tetravalent (Zr, C, and Si), and pentava- 
lent (V and Nb) impurities in LiFeP04 and possible de- 
fect complexes between the extrinsic impurities and na- 
tive defects. The effects of the impurities on the ionic and 
electronic conduction will be discussed and comparisons 
with available experimental data will be made where ap- 
propriate. On the basis of our studies, we discuss spe- 
cific strategies for enhancing the electrical conductivity 
in LiFeP0 4 . 



II. METHODOLOGY 

Our first-principles calculations of native defects and 
extrinsic impurities (hereafter commonly referred to as 
"defects" ) were based on density-functional theory within 
the GGA+U framework) 1 ^— which is an extension of 
the generalized-gradient approximation^ and the pro- 
jector augmented wave method ; 23 ' 2 as implemented in 
the VASP codei 2 ^— The U value for Fe was taken from 
Zhou et ali2&. We used an orthorhombic (1x2x2) super- 
cell, corresponding to 112 atoms per cell, to model the 
defects in LiFeP04. 

Different defects in LiFeP0 4 are characterized by their 
formation energies. The formation energy (E^) of a de- 
fect is a crucial factor in determining its concentration. 
In thermal equilibrium, the concentration of the defect 
X at temperature T can be obtained via the relation^i^ 

c(X) = Ar sites jV co ^ g exp[-£'(X)/fc s T], (1) 

where AT sites is the number of high-symmetry sites in the 
lattice per unit volume on which the defect can be incor- 
porated, and iVconfig is the number of equivalent config- 
urations (per site). It emerges from Eq. ([TJ that defects 
with lower formation energies will occur in higher con- 
centrations. 

The formation energy of a defect X in charge state q 
is defined as^i& 

^(X 9 ) = E tot (X«) - fit* (bulk) -J2 n ^ 

i 

+q{E v + AV + Me ), (2) 

where £? to t(X 9 ) and i? to t(bulk) are the total energies of 
a supercell containing X and of a supercell of the per- 
fect bulk material; //, is the atomic chemical potential of 
species i (and is referenced to the standard state) , and m 
denotes the number of atoms of species i that have been 
added (rii>0) or removed (rij<0) to form the defect. /i e 
is the electron chemical potential (the Fermi level) refer- 
enced to the VBM in the bulk (E v ). AV is the "potential 
alignment" term, i.e., the shift in the band positions due 
to the presence of the charged defect and the neutraliz- 
ing background, obtained by aligning the average elec- 
trostatic potential in regions far away from the defect to 
the bulk valued 
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The atomic chemical potentials fXi are variables and 
can be chosen to represent experimental conditions, 
and are subject to various thermodynamic constraints. 7 
These constraints ensure that LiFePC>4 is thermodynam- 
ically stable. The calculated phase diagrams of the qua- 
ternary Li-Fe-P-02 system at K reported by Ong et 
ali^E show that the compound is stable over a range of 
the oxygen chemical potential (po 2 ) values, from —3.03 
to —8.25 eV. These two values are the upper and lower 
limits in the range of fj,o 2 values considered in our work. 
Lower /io 2 values represent the so-called "more reducing 
environments," which are usually associated with higher 
temperatures and lower oxygen partial pressures and/or 
the presence of oxygen reducing agents; whereas higher 
/io 2 values represent "less reducing environments."— For 
each /xo 2 , the chemical potential of Li, Fe, and P can 
be defined through a set of constraints relevant at that 
/io 2 value. For example, at no 2 =— 4.59 eV, the ther- 
modynamically allowed range of /LtLi and /ip e values is 
the area defined by secondary phases Fe2C>3, Fe 3 (PC>4)2, 
Fe2P2C>7, Li4P2C>7, and Li3PC>4 in the chemical-potential 
diagram as depicted in Fig. 1 of ReL 7 -; the remaining 
variable fip is determined via the condition that ensures 
the stability of LiFeP04f£ Our investigation of thermo- 
dynamically allowed /Un, HF e , /ip, and [io 2 values thus in- 
volves all possible Li-Fe-P-02 phases reported in ReL22. 
Environments under which the system is close to form- 
ing Li-containing (Fe-containing) secondary phases are 
referred to as Li-excess (Li-deficient) iL22 

For the extrinsic impurities, the lower limit on /i.; is 
minus infinity and the upper limit is zero (with respect 
to the energy of the elemental bulk phase) J£ Although 
stronger bounds on the impurity chemical potentials can 
be estimated based on other solubility-limiting phases 
formed between the impurities and the host constituents, 
there is no way to know their exact values. In the fol- 
lowing presentation, the atomic chemical potentials of 
the impurities were chosen somewhat arbitrarily within 
the range given by the above mentioned upper and lower 
limits, provided that the calculated formation energies 
at /Ltg nt are positive. These choices, however, in no way 
affect the physics of what we are presenting since we are 
interested only in the relative formation energies of the 
impurities at a given set of Li, Fe, P, and O2 chemical 
potentials. In other words, we are trying to answer the 
following question: Once a certain impurity is incorpo- 
rated into the material, what is its lattice site preference 
under thermodynamic equilibrium conditions and what 
is its possible effect on the material's properties? Further 
details of the calculations, our methodology (and its lim- 
itations), and the results for native defects in LiFeP04 
can be found in Refi 



III. RESULTS AND DISCUSSION 

Impurities in LiFeP04 can be substitutional at the Li, 
Fe, and P sites, or may be at interstitial sites. For each 



substitutional or interstitial impurity, calculations were 
carried out in several possible charge states. We paid 
attention to positively and negatively charged impuri- 
ties that have the lowest formation energy at /j,™*, and 
hence the highest concentration, since they can act as 
acceptor-like (donor-like) dopants that shift the Fermi 
level toward the VBM (CBM). Note that /4 nt is depen- 
dent on the chemical potentials. For example, we find 
that /4 nt =0.58-0.59 eV for ^o 2 =-3.03 eV, 0.72-1.06 
eV for ^o 2 =-4.59 eV, and 1.98-2.00 eV for fi 0a =-8.21 
eVi£ We also considered possible complexes consisting 
of extrinsic impurities and native point defects. In the 
following, we present only results for those defect config- 
urations that have low formation energies and hence are 
potentially relevant to the material's properties. Note 
that the defect formation energies presented in all the 
figures were evaluated at /xo 2 =— 3.03 eV (corresponding 
to, e.g., a temperature of 800°C and an oxygen par- 
tial pressure of 10~ 3 atm, Refi^S), and equilibrium be- 
tween LiFcP04, Fe 2 03, and Fe7(P04)6 was assumed (i.e., 
Li-deficient environment)^ These conditions give rise to 
/j,Li=— 3.41 eV, fiFc=— 3.35 eV, and [ip=— 6.03 eV with 
the formation enthalpies of the Li-Fe-P-02 phases taken 
from Refi22. Under these conditions, p + and VjJ are the 
dominant native point defects and the Fermi level is at 
/4 nt =0.59 eVr as illustrated in Fig. Q] Of course the 
energy landscape may be different for a different set of 
Li, Fe, P, and O2 chemical potentials. Therefore, we will 
extend our discussions and conclusions to all thermody- 
namically allowed chemical potentials mentioned in the 
previous section. 



A. Monovalent impurities 

Figures HJa) and [DJb) show the formation energies of 
Na- and K-related substitutional impurities at the Li and 
Fe sites, evaluated at ^o 2 =^3.03 eV and its associated 
MLi, A*Fc, and ^p values as noted earlier. The slope in the 
formation energy plots indicates the charge state. Posi- 
tive slope indicates that the defect is positively charged, 
negative slope indicates the defect is negatively charged; 
cf. Eq ^j. We find that Na^; and Nap e are elementary 
point defects; whereas Na^; is, in fact, a complex of Na^ 
and small hole polaron p + that is stabilized at the neigh- 
boring Fe site, and Na Fo a complex of Nap e and p + . At 
the Fe site where the hole polaron is located, the aver- 
age Fe— O bond length is 2.06 A (compared to 2.18 A of 
the other Fe-0 bonds) and the calculated magnetic mo- 
ment is 4.28 fiB (compared to 3.76 at other Fe sites); 
these numbers are in excellent agreement with those ob- 
tained for an isolated hole polaron in LiFeP04. 7 For all 
possible Li, Fe, P, and O2 chemical potentials, Na^j has 
the lowest formation energy at //J 1 *, as seen in Fig. UK a) 
for fi . 2 —— 3.03 eV, suggesting that Na is likely to exist 
as neutral Na^ under thermodynamic equilibrium and 
therefore not effective in enhancing either p + or VjJ . The 
formation energy difference between Na^; and Nap e is, 
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FIG. 2: Formation energies of (a) Na- and (b) K-related im- 
purities at /io 2 — — 3.03 eV, plotted as a function of Fermi level 
with respect to the VBM. 



however, small (about 0.25 eV at most, as measured at 

Strictly speaking, Naj^ is not the positive charge state 
of substitutional Na at the Li site, and Na Fc is not the 
neutral charge state of Na at the Fe site. In other words, 
NaLi, as a single point defect, is only stable as Na^i, and 
NaFe as Na Fe , since locally stable configurations of charge 
states other than Na°; and Na Fc cannot be stabilized. 
If we try, e.g., to create Najj, it decays to a situation 
where the positive charge is not associated with the point 
defect but corresponds to a small hole polaron stabilized 
at the neighboring Fe site, as presented above. Note that, 
instead of writing down, e.g., Na^ as NaL;-p + , we keep 
the notation Na^; for simplicity, but it should be regarded 
as a nominal notation referring to a defect complex that 
consists of the stable/elementary Na^ and a small hole 
polaron. 

For K-related impurities, K°j and K Fe are elementary 
defects; whereas K^; is a complex of and p + , and 
K Fe a complex of K Fe and p + , all similar to Na- related 
impurities. However, we find that, for /io 2 =— 3.25 eV 
and Li-excess or for —3.25 eV</zo 2 <— 3.03 eV, K Fc and 
have the lowest energies for the entire range of the 



Fermi-level values, similar to the energy landscape pre- 
sented in Fig. [U(b); whereas (or K Fc , depending on 
the specific set of the chemical potentials) has the lowest 
energy for /io 2 =™3.25 eV and Li-deficiency or for —8.21 
eV</xo 2 <~ 3.25 eV. The fact that K Fe can have lower 
formation energy than any other K-related impurity at 
/i' e nt suggests that this defect can shift the Fermi level 
toward the VBM and lower the formation energy of pos- 
itively charged native defects. As a result, there will be 
an increase in the concentration of p + and hence in the 
electronic conductivity. It is, however, noted that the for- 
mation energy difference between K Fc and K° ; (or K Fo ) 
is small (about 0.15 eV at most). The difference between 
the general results for Na and those for K, albeit small, 
can be attributed to the difference in their ionic radii. 
Note that K can be incorporated at the Fe site only for 
a small range of the oxygen chemical potential (—3.25 
eV</xo 2 <— 3.03 eV) under thermodynamic equilibrium, 
and the Fermi-level shift is relatively small (about 0.18 
eV at most). The effects of K, therefore, may not be very 
significant. 

Similar calculations were also carried out for Cu- and 
Ag-related impurities. We find that Cu^ has the lowest 
formation energy for /io 2 <— 4.59 eV, whereas Cu Fo has 
the lowest energy for ^o 2 >— 4.59 eV. For Ag-related im- 
purities, Ag°; has the lowest formation energy at /i™* for 
all thermodynamically allowed Li, Fe, P, and O2 chemi- 
cal potentials. Therefore, Cu and Ag are not likely to be 
effective in shifting the Fermi level, similar to Na. 

Experimentally, it has been suggested that Na- and 
K-doped LiFeP04/C samples have an improved elec- 
trochemical performance compared to the undoped 
LiFeP0 4 /C^i^ For example, Na doping was re- 
ported to slightly enhance the electrical conductiv- 
ity: 1.9xl0~ 2 and 0.55xl0~ 2 S.cm -1 for, respectively, 
Li .97Nao.o3FeP04/C and undoped LiFePO^/C at room 
temperature^ As discussed above, Na at the Li site is 
electrically inactive and thus cannot enhance the con- 
centration of p + or V{[ and hence the conductivity. Fur- 
thermore, Na is less likely to be incorporated at the Fe 
site. One, however, cannot exclude the possibility of also 
having a significant concentration of Na at the Fe sites, 
especially if the samples are prepared under conditions 
far from equilibrium. Given the relatively small forma- 
tion energy difference between the Li and Fe sites, fur- 
ther experimental studies should be carried out to clarify 
the lattice site preference of the monovalent impurities in 
LiFeP04 and its dependence on the synthesis conditions. 



B. Divalent impurities 

Figure [3] shows the formation energies of Mg-related 
substitutional impurities at the Li and Fe sites, evalu- 
ated at ^o 2 =-3.03 eV. We find that Mg Fc and Mg^ are 
elementary point defects, whereas Mgff is a complex of 
MgJ; and p + . Among all possible Mg-related impurities, 
Mg Fc has the lowest formation energy at /Zg nt for all pos- 
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FIG. 3: Formation energies of Mg-related impurities at 
fJ.o 2 ~— 3-03 eV, plotted as a function of Fermi level with re- 
spect to the VBM. 
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FIG. 4: Formation energies of Al-related impurities at 
fio 2 — — 3.03 eV, plotted as a function of Fermi level with re- 
spect to the VBM. 



sible Li, Fe, P, and O2 chemical potentials. Note that 
the formation energy at /z™* of Mg^; is higher than that 
of Mg Fc by 0.13—0.36 eV, depending on the specific set 
of the chemical potentials. In fact, the energy of Mg Fc 
is lowest for almost the entire range of the Fermi-level 
values except very near the VBM, as seen in Fig. [3] for 
Ho 2 =— 3.03 eV. Thus Mg is likely to be electrically inac- 
tive. Zn gives similar results (not shown in the figure), 
i.e., Zn is likely to exist as neutral Zn Fc and does not 
shift /i e away from /i™* . 

Experimentally, Roberts et al.— reported that there 
was no evidence of Mg on the Li site in samples prepared 
with the nominal stoichiometry Lii-^Mgo.szFePO^ 
which is consistent with our conclusion that Mg is likely 
to exist as Mg Fc in LiFcPC>4. Note that these results are 
in contrast to those reported by other research groups 
who suggested that Mg occurs on the Li site^ ' 34 ' 35 Since 
Mg Fc can enhance neither the concentration of p + nor 
that of VjJ, the reported enhancement in the electrical 
conductivity of Mg-doped LiFeP0 4 (Refs^^) is dif- 
ficult to understand as arising directly from Mg doping 
at the Fe site. However, given the relatively small for- 
mation energy difference at /i' e nt between the Li and Fe 
sites (0.13—0.36 eV), a significant concentration of Mg 
at the Li sites might still be possible, especially if the 
samples are prepared under conditions far from equilib- 
rium. Regarding Zn, Bilecka et ali^I reported a con- 
ductivity of 10.42xl0 -8 S.cm -1 at room temperature in 
LiZn x Fei_a;P04 samples, compared to 7.5xl0~ 8 S.cm -1 
of undoped LiFeP04. Since Zn Fc is not effective in shift- 
ing the Fermi level, this small conductivity increase may 
be due to other effects. 



C. Trivalent impurities 

Figure 2] shows the formation energies of Al-related 
substitutional impurities at the Fe and P sites, evaluated 
at no 2 — ~3.03 eV. We find that Al Fc is an elementary 
point defect, whereas Alp^t is a complex of Al Fe and p + . 
At the P site, Alp~ is found to have the lowest energy, 
which is not unexpected, given the valence of Al (+3) and 
P (+5). In this configuration, Al replaces P and forms a 
slightly distorted AIO4 unit with the average Al— O bond 
length of 1.77 A, compared to 1.56 A of the P— O bonds. 
Among the Al-related impurities, we find that Al Fe has 
the lowest formation energy at /z^ for all thermodynam- 
ically allowed Li, Fe, P, O2 chemical potentials. The for- 
mation energy of the lowest energy configuration of Alu 
is found to be higher than that of Alp e by 0.65—0.89 eV, 
as measured at ^™ t . 

With Al existing as the positively charged impurity 
Alp e , more V£[ and Li Fe are created to maintain charge 
neutrality, thus the Fermi level is shifted away from /i' e nt 
and toward the CBM. As a result of this shift, the for- 
mation energy (concentration) of p + is increased (de- 
creased). We also investigated possible defect complexes 
between Al Fe and the negatively charged native point 
defects, i.e., Alp^-V^ and Al Fc -Li Fc . We find that the 
binding energy of these complexes are, respectively, 0.38 
and 0.25 eV with respect to their isolated constituents. 
For any given set of the atomic chemical potentials, Al Fc - 
Lip e has a slightly higher (0.04—0.23 eV) formation en- 
ergy than Alp c -T^ L T. The distance between the two defects 
in Alp c -Lip c is 3.28 A, whereas it is 3.96 A in Al Fe -Li Fe . 
Our results suggest that, although Al Fc enhances both 
and Lip c , the concentration enhancement is larger 
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FIG. 5: Formation energies of Zr-related impurities at 
/io 2 — — 3.03 eV, plotted as a function of Fermi level with re- 
spect to the VBM. 

for the vacancies than for the antisites. 

Experimentally, Amin et aL 38 ' 39 reported that Al- 
doped LiFeP04 has a higher ionic conductivity but lower 
electronic conductivity compared to undoped LiFeP04. 
This is exactly consistent with our results that Al Fe en- 
hances and reduces p + . Regarding the lattice site 
preference, although Amin et al. suggested that Al oc- 
curs on the Fe site, which is in agreement with our re- 
sults, other research groups reported that the dopant is 
predominantly on the Li site4^i2£ Further experimental 
studies, therefore, should be carried out to clarify this 
situation. 



D. Tetravalent impurities 

Figure [5] shows the formation energies of Zr-related 
substitutional impurities at the Fe and P sites, evaluated 
at no 2 =— 3.03 eV. Potentially relevant defect configura- 
tions are Zr F +, Zr F + (a complex of Zr F + and p + ), and 
Zrp. In Zrp, Zr replaces P and forms a ZrC>4 unit with 
the average Zr— O bond length of 1.98 A, compared to 
1.56 A of the P-0 bonds. We find that Zr 2 + has the 
lowest formation energy at ^ nt for all possible Li, Fe, P, 
O2 chemical potentials, suggesting that Zr has a donor- 
like effect, i.e., shifting the Fermi level toward CBM and 
hence enhancing the negatively charged native defects. 
The formation energy of the lowest energy configuration 
of ZrLi at //{, nt is found to be higher than that of Zr Fe by 
1.17-1.41 eV. 

We also considered complexes between the elementary 
Zr F + and the negatively charged native point defects. 
Zr F +-2V L T is found to have a binding energy of 0.78 eV 
with respect to its isolated constituents; the distances 



between Zr F + and the vacancies are 3.26 and 3.62 A. 
Zr F +-2Li Fe has a binding energy of 0.65 eV; the distances 
between Zr F + and the antisites are 4.05 and 4.13 A. Un- 
der Li-deficient environments, the formation energy of 
Zrp^-21/jT is found to be lower than that of Zrp^"-2Li Fe 
by about 0.20—0.34 eV as seen in Fig. [5] whereas un- 
der Li-excess environments, the energy of Zr F ^~-2V^T is 
comparable to that of Zr F +-2Li Fc . Our results therefore 
suggest that Zr F + is more likely to be associated with 
Vj^ than with Li Fe , thus enhancing the vacancies and 
the ionic conductivity. 

The donor-like effect of Zr in LiFeP04 could have con- 
tributed to the enhancement in the total conductivity 
as observed in experiment.— However, in order to clearly 
see the effect, one should carefully look at the ionic and 
electronic conductivities in undoped and doped samples 
(preferably without carbon coating). Note that, an en- 
hancement in the intrinsic ionic conductivity usually re- 
sults in a reduction in the intrinsic electronic conductiv- 
ity, as discussed in the Introduction. Regarding the lat- 
tice site preference, several experimental works suggested 
that Zr occupies the Li site^^- which is in contrast to 
our results showing that the Fe site is energetically more 
favorable under thermodynamic equilibrium conditions. 

Figures EJa) and EUb) show the formation energies of 
C- and Si-rclatcd substitutional impurities at the P site, 
evaluated at fio 2 = ~ 3-03 eV. Cp is an elementary defect, 
whereas C P is a complex of Cp and p + , Cp a complex 
of Cp and two p + . We also considered Cp-Vo + , a posi- 
tively charged complex consisting of Cp and an oxygen 
vacancy (Vq + ). This complex has a binding energy of 
4.44 eV with respect to its isolated constituents. It can 
also be regarded as replacing a (P04) 3- tetrahedron with 
a trigonal planar (COa) 2 ". In the (COs) 2- unit, the av- 
erage C— O bond length is 1.31 A, compared to 1.56 A 
of the P— O bonds. C thus exhibits a +4 valence at the 
P site. At the Li and Fe sites, we however find that C 
exhibits a +2 valence. As a result, C^; and C Fe are the 
elementary defects, whereas their higher charge states are 
defect complexes consisting of (or C Fe ) andp + . 

For -3.25 eV</i O2 <-3.03 eV, Cp~ has the lowest for- 
mation energy at /i' e nt . Given its positive effective charge, 
this defect would exhibit an donor-like effect, enhancing 
the concentration of negatively charged native defects. 
Moreover, for each Cp created, there are two p + ; thus, 
with Cp being the lowest energy configuration, C can 
enhance both p + and (see more below). For —8.21 
eV<no 2 <~ 3.25 eV, the formation energy of C F -V^ + is 
lower than that of Cp , and the complex has the lowest en- 
ergy at /Zg nt . Since Cp-V^" 1 " exhibits a donor- like effect, it 
can enhance VjJ. In fact, we have investigated a complex 
of Cp-Vo and Vj^, hereafter denoted as Cp-Vo "^lp 
and find that it has a binding energy of 0.63 eV with re- 
spect to its isolated constituents. The formation energy 
of this complex [see Fig. [6ja)] is lower than that of C F - 
(not shown in the figure), a neutral complex of Cp 
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FIG. 6: Formation energies of (a) C- and (b) Si-related impu- 
rities at (J,o 2 =— 3-03 eV, plotted as a function of Fermi level 
with respect to the VBM. 



and VjT, by 0.25—3.79 eV. Our results suggest that, in 
the presence of V^J, Cp-V^ -1- is energetically more favor- 
able than Cp for all possible Li, Fe, P, and O2 chemical 
potentials. In other words, if incorporated into the bulk 
of LiFePCU, C is likely to exist as Cp-V^-V^, and the 
presence of C thus enhances V{~± only. The formation 
energy of C at the Li and Fe sites is much higher than 
at the P site, with the energy difference at ^t™* between 
Cp and Cfc is 4.80—6.24 eV. Note that, here we assume 
C can go into the bulk with a concentration higher than 
that of the native defects. The penetration of C, however, 
may be impeded by reactions at the surfaced Since C 
is present in precursors used in the synthesis of LiFeP04 
and also because LiFeP04 is usually coated with C, it is 
worth carrying out further studies to see if there is C in 
the interior of the material. 

Regarding Si-related substitutional impurities at the P 
site, Sip is found to be an elementary defect, whereas Si P 
is a complex of Sip and p + , and Sip a complex of Sip 
and two p + . Sip-V^ + has a binding energy of 2.35 eV 
with respect to Sip and Vq + . In this complex, Si forms 
a SiO.4 unit that shares one O atom with the neighbor- 
ing PO4 unit, i.e., Si03-0-P03 with the Si— P distance 



being 2.90 A, compared to 3.71 A of the P— P distance. 
Energetically, we find that Sip has a lower formation en- 
ergy than Si p -Vo + except for -8.21 cV<^o 2 <-7.59 eV. 
Sip is found to have the lowest formation energy at fi™ 
for —4.59 eV<^o 2 <— 3.03 eV, as seen in Fig. [SJ^b) for 
fJ.o 2 — ~ 3-03 eV. Although Sip, is neutral and thus not ef- 
fective in shifting the Fermi level, each Si P contains one 
p + . The formation of Si P , therefore, may help introduce 
more p + into the material. For —8.21 eV</j,o 2 <— 4.59 
eV, Sip is found to have the lowest formation energy. 
In this case, Si exhibits an acceptor-like effect and is ex- 
pected to enhance the concentration of p + . For any given 
set of Li, Fe, P, and O2 chemical potentials, neutral com- 
plexes such as those consisting of and Sip or Sip-VQ + 
all have higher formation energies than Si p . Our results 
thus suggest that Si incorporated at the P site enhances 
small hole polarons. At the Li and Fe sites, Si exhibits 
a donor-like effect since Si is supervalent to Li and Fe; 
the elementary, lowest energy configurations at /Xg nt are 
Sij^ and Si F ^ , where the formation energy of Si at the Li 
site is higher than that at the Fe site by 0.43—0.65 eV, 
depending on the specific set of the chemical potentials. 
The formation energy difference at /j,™' between Sip and 
Si Fc is 0.33-3.56 eV 

Amin et al4^ reported that doping LiFeP04 with Si 
results in an increase of the ionic conductivity and a de- 
crease of the electronic conductivity, suggesting a donor- 
like effect. As discussed above, this effect should be as- 
sociated with Sip^. Our results, however, suggest that 
Si F ^" is not the lowest defect configuration under thermo- 
dynamic equilibrium conditions. Although the cause of 
this discrepancy between experiment and our results on 
the lattice site preference of Si is still not clear, we suspect 
that doping at the P site may be kinetically hindered, 
and thus a thermodynamic equilibrium approach is not 
applicable*^ This may be the case for all dopants that in- 
volve P-site doping in L1FCPO4. It should be noted that 
LiFeP04 can be regarded as consisting of Li + , Fe 2+ , and 
(P04) 3 ~ with ionic bonding between the units. Within 
the (P04) 3 ~ unit, however, the P— O bonds are highly 
covalent. 7 If doping occurs via diffusion, P-site doping 
would involve diffusion of P 5+ and the dopant which is 
expected to be difficult because it necessarily involves 
breaking/forming of chemical bonds between O and P 
(and the dopant). Thus P-site doping might also be af- 
fected by the P- and dopant-containing precursors used 
in the synthesis. 



E. Pentavalent impurities 

Figures EJa) and [7](b) show the formation energies of 
V- and Nb-related substitutional impurities at the Li, 
Fe, and P sites, evaluated at ^o 2 — — 3.03 eV. We find 
that V P is an elementary point defect, suggesting that 
V (vanadium, not to be confused with a vacancy, V) at 
the P site has a +5 valence. Vp is a complex of V P and 
p + , whereas Vp is a complex of V P and a small electron 
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polaron (p~)£ At the Li and Fe sites, however, V exhibits 
a +3 valence, with V L + and Vp e being the elementary 
defects. Attempts to create higher charge states result 
in complexes consisting of these defects and p + , whereas 
lower charge states are complexes containing jT . We 
have also investigated neutral complexes such as Vp e - 
V{^, a complex of Vp~ e and a Li vacancy, and 2Vp c -V F 2 c ~, 
a complex of two Vp e and an Fe vacancy. Energetically, 
we find that V P has the lowest formation energy at /zj, nt 
for all possible Li, Fe, P, and O2 chemical potentials. 
With this neutral configuration, V would be ineffective 
in shifting the Fermi level. 

V has positive effective charges at the Li and Fe sites 
and thus exhibits a donor-like effect. The formation en- 
ergy at /Xg nt of VFe in its lowest energy configuration is 
higher than that of V P by 0.50-2.14 eV for almost all 
possible chemical potentials, and by 0.10—0.20 eV for a 
small region of the chemical potentials associated with 
fio 2 around —4.59 eV. On the other hand, the forma- 
tion energy of Vf is lower than that of Vu by 0.76—0.97 



eV for almost all possible chemical potentials, except for 
a small region associated with ^o 2 <— 4.59 eV where the 
two defects have comparable formation energies. Neutral 
complexes such as Vp^-V^J and 2V Fe -Vp J e ~ have formation 
energies higher than V P by 0.36-2.32 and 4.51-9.09 eV, 
respectively. Again, these results suggest that V is en- 
ergetically more favorable at the P site under thermody- 
namic equilibrium conditions. 

Experimentally, Hong et al42 reported that their at- 
tempt to incorporate V into LiFePC>4 at the Fe site 
resulted in having V at the P site instead. However, 
Omenya et al.— later reported that the substitution at 
the P site could not be reproduced, but at least 10 mol% 
of the Fe sites were occupied by V 3+ . Zhang et al.— also 
reported a valence between +3 and +4 for V in LiFeP04. 
Note that, according our results dicussed above, the +3 
valence of V is associated with substitution at the Fe 
site that exhibits a donor-like effect. Thus, for V-doped 
LiFePC>4 samples reported in Refs . 43 ' 44 , V may indeed 
have been incorporated at the Fe site, although this site 
is found to be energetically less favorable in our calcu- 
lations. The discrepancy between experiment and our 
results on the lattice site preference of V could be that 
P-site doping is thermodynamically favorable but kinet- 
ically hindered, as discussed earlier. 

Regarding Nb, we find that Nb P is an elementary point 
defect, and the dopant thus has a +5 valence at the P 
site. Nb P is a complex of Nb P and p + , whereas Nb P is a 
complex of Nb p and p~ . Nb also exhibits the +5 valence 
at the Fe site, with Nbp^~ being the elementary defect. 
At the Li site, we however find that Nb has a +4 valence 
and the elementary configuration is Nb 3 ^ 1 " ■ Higher charge 
states such as Nb 4 ^ is a complex of Nb 3 ^ and p + . En- 
ergetically, Nb p has the lowest formation energy at /j,™* 
for almost all possible chemical potentials, except for a 
very small region associated with —4.59 eV<no 2 <—3.89 
eV of the (/iLi, MFe, Mo 2 ) polyhedron that defines the sta- 
ble range of Li, Fe, P, and O2 chemical potentials (Ref. 7 ) 
where Nb F ^~ has the lowest formation energy. For exam- 
ple, in the ^o 2 =~4.59 eV plane, the small region near 
Points B, C, and D in the chemical-potential diagram as 
depicted in Fig. 1 of Reffi is where Nbp^" is most favorable 
energetically. Nb P cannot enhance the concentration of 
either or p + , but Nbp^ is donor- like and thus can 
enhance V^T. We also find that Nbp^-SV"^, a complex 
of Nbp^" and three Li vacancies, has a lower formation 
energy than Nb P in that small region of the (/iLi, MFc, 
Mo 2 ) polyhedron. Our results thus suggest that Nb can 
be incorporated at the Fe or P site, depending on the 
synthesis conditions. The reported enhancement in the 
total conductivity of Nb-doped LiFeP04 (Refi 4 -) could 
have been partly due to the enhancement of V L ~[ caused 
by Nbp^". In addition, the incorporation of Nb at the P 
site may be kinetically hindered as discussed above in the 
case of Si and V. 

Overall, we find that extrinsic impurities in LiFeP04 
each have one stable charge state (i.e., elementary defect 
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configuration) at a given lattice site. Any attempt to cre- 
ate higher (lower) charge states will result in complexes 
consisting of the elementary configuration and small hole 
(electron) polarons, which is similar to what has been 
observed for native defects in LiFePC>4 as reported in 
Ref.— . Our studies help identify specific impurities that 
are likely to be effective (or ineffective) for enhancing 
the concentration of p + or VT, . The results also suggest 
further experiments to investigate the lattice site pref- 
erence (and its dependence on the synthesis conditions) 
and changes in the ionic and electronic conductivities. 
For practical applications, since doping alone cannot en- 
hance both p + and VjJ as needed for high ionic and elec- 
tronic conductivities, other methods such as carbon coat- 
ing (Refs i 45 ^ 46 ) and/or thermal treatment (Ref4£) may 
still be needed. As discussed by Julien et al4S, carbon 
coating can help repair structural damage at the surface 
and ensure the electric contact between LiFePC>4 parti- 
cles. Also, with whatever method is employed to enhance 
the conductivity, one should first be able to control the 
experimental conditions during synthesis to reduce Fe^; 
and enhance p + and X 

Note that computational studies of impurities in 
LiFePC>4 have also been carried out by Islam et al . 48 ' 49 . 
Based on calculations using interatomic potentials, they 
reported that, among other impurities, Na is energeti- 
cally more favorable at the Li site, whereas Mg, Al, Zr, 
and Nb are more favorable at the Fe site. These results 
are generally in qualitative agreement with our results 
for these substitutional impurities at the Li and Fe sites. 
However, they did not report results for the impurities 
at the P site, and also did not explore all possible Li, 
Fe, P, and O2 chemical potentials. As discussed earlier, 
the lattice site preference can be sensitive to the chemical 
potentials (e.g., in the case of K and Cu), and the Li and 
Fe sites may not be energetically most favorable under 
thermodynamic equilibrium. 



namically allowed atomic chemical potentials, Na, Cu, 
Ag, Mg, and Zn are likely to exist as neutral defects 
in LiFeP04 and thus do not enhance the concentration 
of either small hole polarons (p + ) or lithium vacancies 
(VJJ), i.e. electronic or ionic conduction. K may be in- 
corporated at the Fe site under thermodynamic equilib- 
rium for a small range of the chemical potentials, where 
it exhibits an acceptor-like effect and thus can enhance 
the concentration of p + . Al and Zr, on the other hand, 
exhibit a donor-like effect and are thus effective in en- 
hancing the concentration of V{^. C can also exhibit a 
donor-like effect if incorporated into the material. We 
find that Si is energetically more favorable at the P site 
where it exhibits an acceptor-like effect, which is in con- 
trast to experiment suggesting substitution at the Fe site 
that exhibits a donor-like effect. V is found to be more 
favorable at the P site and thus electrically inactive, in 
contrast to experiment suggesting substitution at the Fe 
site. Nb is also found to be electrically inactive at the 
P site, except for a small range of the atomic chemical 
potentials where it is energetically more favorable at the 
Fe site and exhibits a donor-like effect. We suggest that 
doping LiFeP04 with Si, V, or Nb at the P site is thermo- 
dynamically favorable but kinetically hindered, and that 
Si, V, and Nb may have been incorporated at the Fe site 
where they exhibit a donor-like effect as reported in ex- 
periment. Indeed, the valence reported for V in V-doped 
LiFeP04 by experiment is consistent with substitution 
at the Fe site in our calculations. Our studies, therefore, 
can serve as guidelines for experiment on which lattice 
site the impurities are most likely to be incorporated un- 
der thermodynamic equilibrium conditions, and how the 
impurities may affect the ionic and electronic conduction. 
To achieve a high electrical conductivity, however, one 
should combine the incorporation of electrically active 
impurities with other methods such as defect-controlled 
synthesis, carbon coating, and/or thermal treatment. 



IV. SUMMARY 

We have carried out first-principles studies of the ef- 
fects of various extrinsic impurities on the ionic and elec- 
tronic conduction in LiFeP04. We find that the forma- 
tion energy and lattice site preference of the impurities 
depend on Li, Fe, P, and O2 chemical potentials which 
represent the synthesis conditions. For all thermody- 
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